o 



(N 



Effective Hamiltonian Crystal Field for 
Magnetic Interactions in Polynuclear 
Transition Metal Complexes. Sequential 
Derivation and Exemplary Numerical 

Estimates. 

M-t' 
g 

(D ! Andrei L. Tchougreeff 

°- ■ March 26, 2013 

o 

>>; 

Moscow Center for Continuous Mathematical Education, Bol. Vlasevskiy 
ft ; per., 14, Moscow 119002, Russia 

(N 

VjQ ■ Abstract 
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I By this we extend our work of the year 1992 devoted to calculating 

the intrashell excitations in the d-shells of coordination compounds of 
the first transition metal row, which resulted in the Effective Hamil- 
tonian Crystal Field (EHCF) method, to their polynuclear analogs in 
order to assure the description of several open d-shells and of mag- 
netic interactions of the effective spins residing in these shells. This 
is a challenging task since it requires improving the precision of ca. 
Ijjj ' 1000 cm -1 (that of describing the excitation energies of the single d- 

■ shells by the already well successful EHCF method) to the that of ca. 

10 -T- 100 cm -1 characteristic for the energies required to reorient the 
spins i.e. eventually by two orders of magnitude. This is performed 
within the same paradigm as used for the EHCF method: the con- 
certed usage of the McWeeny's group-function approximation and the 
Lowdin partition technique. These are applied to develop the effective 
description of the d-system of the polynuclear complexes, composed of 
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several d-shells, including the working formulae for the exchange pa- 
rameters between the <i-shells belonging to different transition metal 
ions. These formulae are implemented in the package MagAixTic and 
tested against a series of binuclear complexes of trivalent Cr and Fe 
cations featuring ^-oxygen superexchange paths in order to confirm 
the reproducibility of the trends in the series of values of exchange 
parameters for the compounds differing by the chemical substitutions 
and other details of composition and structure. The results of calcula- 
tions are in a reasonable agreement with available experimental data 
and other theoretical methods. 

1 Introduction 

Magnetic properties of transition metal complexes (TMC) derive from the 
classical work reducing them to those of the respective isolated d- shells. [T] 
Those of so called polynuclear TMC (PTMC - those containig two or more 
transition metal ions - TMI's - with open d-shells) in the low-energy range are 
usually as well classically mapped to the effective exchange (spin) Dirac-Van 
Vleck-Heisenberg Hamiltonian: [2j 

-^DvVH = 2_j JifiiSj, (1) 
i<j 

controlled by the effective exchange parameters J;.,- describing the interac- 
tions of the effective spins Si residing in the d-shells of the corresponding 
TMI's with those in their neighbourhood. Numerous attempts have been 
undertaken to derive reliable numerical values of these parameters suitable 
to describe the lowest-energy range of excitations of the electronic subsystem 
of PTMCs quantum chemically. The complexity of this task (examples of the 
tentative application objects from the recent literaturejUHUEJE] are quite im- 
pressive) prevented, however, from the decisive success. The theoretical con- 
siderations of the effective magnetic parameters present in the literature are 
numerous, but can be classified in three types: (i) ab initio approaches, (ii) 
DFT-based approaches, and (iii) semiempirical/semiquantitative approaches 
predominantly stemming from the solid state physics realm. 

The ab initio methods are potentially capable to solve the problem (see 
below). (71 El HI [TOl [II] However, the systems of real interest in this area may 
contain thousands of atoms of which hundreds are going to be transition 
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metal ones. [31 SI El E] Within the ab initio realm such systems require the 
account of correlations (both dynamical and noindynamical) at a level which 
causes the scaling of the requred computational resources as N 7 where N is 
the number of orbitals involved. It will be hardly possible to apply such an 
approach within a chemical problem setting that of establishing/reproducing 
the trends in a series of similar compounds rather obtaining a unique number 
for a unique system. The dependence of the result on the basis state used 
has been established as well. [9] 

The DFT based methods recently received considerable attention in this 
context. [T2| 113] The success reported for the DFT-based techniques [121 H3] 
heavily relies upon the possibility to obtain the broken symmetry solutions in 
the UHF setting for the respective molecules. Although trivial for binuclear 
complexes, usual guinea pigs for the workers in the field, it may become very 
problematic for the systems where odd-numbered cycles of spins coupled 
through antiferromagnetic exchange interactions can be singled out and thus 
one can expect frustration, so that the spins cannot arrange themseves in 
a single strucutre describable by the Ising model which is a prerequisite for 
applying either the DFT-based technology or similar approach based on usage 
of nonempirical Hamiltonians. [T2| IT3] □ 

These complications made us to reconsider the ingenious idea due to P.W. 
Anderson [161 E] wno suggested a two step procedure eventually yielding 
the reliable estimate of the magnetic exchange parameters without invoking 
a priori unknown magnetic ground states of a complex or a solid. According 
to the Anderson's prescription (see also [181 HI]) i n a line with the original 
implicit picture used by BethejT] formulated explicitly only later. [20] At the 
first step one has to obtain one-electronic wave functions of the magnetic 
ions (in our case these are the states of the TMI's d-shells) immerced in 
the insulating diamagnetic medium (ligands) within some kind the ligand 
field theory [2T| 22J At the second step one has to estimate the interactions 

1 It must be understood, as well, that the above examples of experimentally interesting 
systems [5j [4J EJ are still not accessible for the DFT based methods either. Also one 
has to admit that usage of the DFT methods for this purpose is problematic in a more 
general sence: already reproducing the ground state total spin, one-center d-d spectra 
and their geometry dependence of the mononuclear TMCs represents a problem for them. 
Apparently the precision requirements for estimating the magnetic interactions in the 
polynuclear TMCs are much more tough than for the ground state spin in mononuclear 
TMCs. Thus the DFT-based methods under certain conditions can be tuned to reproduce 
the observable exchange parameters, but the uniform picture of (P)TMCs evades the DFT 
methods. [T4l[T5] 
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between electrons residing in the magnetic orbitals defined at the first stepo 
This way of thinking was pursued in our work|23j where we performed 
the first step: the consistent description of the magnetic states of the iso- 
lated TMIs. It is commonly believed that semi-empirical methods of quan- 
tum chemistry are not of much use for TMCs. The reason of such attitude 
originates from very narrow understanding of the nature of semiempirical 
methods. Of course, if one takes them naively: as a parameterized SCF in 
the valence basis, one cannot count on considerable success. Indeed, such an 
approach is known to be unsuccessful for decades in coping with the TMCs 
ground state total spins and with their geometry dependence not talking 
about the low-energy d-d excitation spectra. The way out had been found|23j 
in the concerted usage of the McWeeny group-function formalism [21] and of 
the Lowdin partition technique[25j as described in details in reviews[26] [27] 
and in monograph. [28J This can be considered as a new concept of semiem- 
pirism in quantum chemistry: taking the form of the wave function on the ba- 
sis of observable electronic groups (like G?-shells, 7r-systems, two-center bonds 
etc.) characteristic for the considered class of molecules. Taking this way 
allowed us previously to reach considerable success in describing the ground 
states and d-d excitation spectra of TMCs[23j [29j E01 EH [32] and of their 
geometry dependence. [HJ EH [331 EH ESI EHJ E7J [38] In the present work we 
concentrate on the derivation of the effective Hamiltonian for the d-system of 
a PTMC which eventually contains necessary effective magnetic terms and by 
this follow the way (iii) of developing a substantiated semiempirical approach 
to estimating the parameters of effective magnetic interactions suitable for 
sizeable systems without addressing in advance their magnetic structure. 
This way would not be possible without works[39j HOI HT1 1351 1121 IM 135] 
where the quantitative empirical relations between the parameters of the 
crystal field felt by the d-shells in series of PTMCs and their respective mag- 
netic exchange parameters has been firmly established. These phenomeno- 
logical works allowed us to hope (and as shown below with good grounds) 
that the estimates of the crystal field parameters by semiempirical procedures 
assembled in the EHCF method|23j apply as well to the magnetic exchange 
parameters in PTMCs. 

The paper is organized as follows. In the next Section we describe the 

2 Calculating system «without magnetism» in a solid state i.e. translationally invariant 
rather a «crystal field» setting would lead to a metallc state with poorly definable Wannier 
states. That is why one has to rely upon as well not well defined cluster models for these 
ions. 
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necessary modification of the EHCF formalism for the case of polynuclear 
TMCs (PTMCs) and specific contributions to the effective Hamiltonian of 
their d-electrons. The details of the derivation are quite cumbersome and 
thus transferred to the Appendices. In further two Sections we describe 
respectively some relevant details of implementation and results of our cal- 
culations as compared to previous ones and the experiment. After that we 
discuss our results and the relation or our approach with others existing in 
this area. 

2 Effective Hamiltonian of Crystal Field as mod- 
ified for magnetic interactions 

The intuitive physical picture we are going to pursue is very simple. The 
low-energy range of excitations of the electronic system of a polynuclear 
transition metal complex (PTMC) is spanned by the states differring by the 
relative orientations of the local spins (electronic in origin) residing in the 
respective d-shells. From the ab initio or DFT viewpoints such a picture is 
an enormous oversimplification since everyone is aware of the fact that all 
the basis electronic states are strongly hybridized. Nevertheless, the relative 
magnitudes of the electron transfer (hopping) matrix elements in compari- 
son to the relevant energies of the one-electron charge transfer states differ 
drastically for different states within the PTMC. Namely, this ratio is by 
an order of magnitude larger for the transfers involving d-AOs than for a 
similar ratio for the transfers involving the AOs on the close ligand atoms. 
This is the physical reason why in most cases one can unequivocally ascribe 
some definite integer number of electrons to the d-shells and interpret the 
low-energy optical spectra of (P)TMCs in terms of the states of such shells 
with fixed integer numbers of electrons i.e. with use of the crystal field 
theory. [2T| [22| H9] These ideas allowed us to make a formal move and to 
identify the parts of systems according to their respective different physical 
regimes: the highly correlated <i-AOs which we term to span the one-electron 
space of the (f-system and the weakly correlated AOs of the «organogenic» 
ligand atoms and the weakly correlated sp-AOs of TMI's which span the 
one-electron space of the /-system. The difference of the physical conditions 
in so defined systems allows us treat them differently as well. The /-system 
being in a low correlated regime can be treated with an acceptable precision 
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by a semiempirical SCF theory (INDO will suffice, see below); the (i-system 
requires in principle kind of correlated treatment. In order to be able to treat 
so defined subsystems separately one has to develop (mutually dependent) 
effective Hamiltonians for the singled out subsystems. That for the d-system 
eventually contains the required «magnetic» terms. This development is 
sketched in the subsequent Subsections. 

2.1 Molecular Hamiltonian breakdown 

The basis of one-electron states for describing a PTMC is defined within a 
general semiempirical scheme and is thus spanned by the valence AOs which 
are considered to be implicitly orthogonal. For a PTMC the valence AO 
basis is formed by the valence AOs of the ligand atoms (largely, but not 
mandatory, organogenic) and by those of the TMIs it contains. These latter 
are the 4s-, 4p-, and 3d- AOs for the atoms of the first transition series. 
On the basis of the ratios of the characteristic intersubsystem hopping and 
Coulomb integrals we separate as explained previously the entire set of the 
AOs into the basis subsets spanning respectively d- and /-systems of a PTMC. 
The corresponding orbital subspaces are thus implicitly orthogonal. 

According to the above classification of the valence AO's the semiempir- 
ical Hamiltonian for a PTMC rewrites as (the explicit form of all operators 
is given in Appendix [Al: 

H = H d + Hi + H mt , , 

H mt = H c + H r + H x 

where the operators Hj (one for <i-electrons of all TMIs in the complex) and 
Hi (that for electrons in the /-system) are formally defined as those involving 
the Fermi creation and annihilation operators[46| referring only to the d- and 
/-systems, respectively, whereas the interaction operator H mt involves mixed 
products of the Fermi creation and annihilation operators referring to the d- 
and /-systems (see Appendix lA.lj) . 

The interaction operator H lTLt is a sum of contributions, which as well 
can be formally defined through the types of the products of the Fermi op- 
erators entering them. Correspondingly, the operators H c and H x are two- 
electron operators which involve the products of four Fermi operators and 
represent, respectively, the Coulomb and exchange interactions between two 
subsystems; the operator H r is the one-electron operator which involves only 
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products of one Fermi operator referring to the d- and one Fermi operator 
referring to the I- AO subspaces (details see in Appendix \A.2\i . 

2.2 Configuration subspaces and Lowdin partitions 

The main complication for a sequential description of magnetic properties 
of PTMCs is the enormous number of electronic configurations to be taken 
into account. Yet a simpler task of describing mononuclear TMCs involved 
quite a number of intermediate («outer space») configurations which were 
addressed with use of the Lowdin partition technique. |25j Its success as well as 
technical details of its implementation in the EHCF method makes us to try 
similar moves as applied to PTMCs. The Lowdin partition can be regarded 
as a specific method of separating variables in a many-particle system such 
that upon applying it, the system turns to be representable by a product 
function whose multipliers describe otherwise entangled states of the parts 
of the whole system. That is to say that performing a Lowdin partition 
is a prerequisite of employing the McWeeny group-function approximation. 
Thus our derivation evolves as a sequence of the Lowdin partitions allowing 
to single out relevant («model») configuration subspaces and to take into 
account the effect of the «outer space» configurations outside the model 
subspaces considered explicitly at each step on the matrix elements of the 
model Hamiptonian. 

2.2.1 Separation of d- and /-electronic variables 

The first move in the EHCF derivation was the separation of the variables in 
the single (f-shell available in a mononuclear TMC and in the ligands. That 
finally allowed us to write the TMCs electronic wave function as an antisym- 
metrized product of the group functions for two distinguishable subsystems 
in the TMC, having physically different conditions, the (f-shell and the rest: 

*TMC = $>d A $ L , (3) 

the form implicitly assumed by the crystal field theory. [H [2U 122^1 The form 
eq. ([3]) is approximate since the operator H r acts and the electrons move 

3 Here A stands for the antisymmetrized product of the electronic wave functions, if 
the Fermi operators are used this notation means that all of them referring to the d-shell 
stand to the left of those referring to the ligand orbitals 
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between the subsystems so that the entire TMC occurs in an electronic state 
where those of the individual subsystems are entangled. In order to recover 
the crystal field theory picture in which the number of electrons in the d- 
system of the TMC is constant for the low-energy part of the electronic 
spectrum the exact Hamiltonian has to be projected to the model subspace 
spanned by the configurations in which the number of electrons in the d- 
shell(s) is fixed. The same move has to be done for the PTMCs. Let the 
overall number rid of electrons in the <i-system be fixed. The numbers of 
electrons in the individual <i-shells satisfy the following conditions: 

m = N-n d , 

where rz^(i) is the number of electrons in the d-shell of the i-th TMI and N 
is the total number of electrons in the PTMC Since the states with different 
distributions of electrons among the c/-shells of the different TMIs enter the 
consideration (see below) it may be technically easier to trace that ni - the 
number of electrons in the /-system - remains constant. As in the derivation 
of the EHCF method|23j the configurations with numbers of d- (or /-) elec- 
trons different from that fixed above (ligand to metal and metal to ligand 
charge transfer states, respectively abbreviated as LMCT and MLCT) have 
to be taken into account to keep track of the effect of the resonance operator 
H r . 

As in the derivation of the EHCF method[23j it is done with use the 
Lowdin partition technique. |25j Let P be the projection operator to the model 
many-electron configuration subspace spanned by the Slater determinants 
with ni electrons in the /-system and remaining rid m ah the <i-shells and let 
Q = 1 — P be its complementary projection operator. The Lowdin partition 
replaces the original Hamiltonian acting in the entire configuration space 
by the effective Hamiltonian H eS (E) which acts in the model configuration 
subspace ImP: 

H cS (E) = PHP + PH„.P, 

but by construction has the same eigenvalues as the original Hamiltonian. In 
eq. ([5]) the resolvent: 

R(E) = (EQ - QHQ)- 1 (6) 

is used. The form and the notation H rr is due to the fact that the operator 
H r , the only term which transfers electrons between the d- and /-systems of 



8 



the PTMC and by this mixes the states in the ImP and ImQ configuration 
subspaces, 

PHP = PH P + PH X P 

Ho = H d + Hi + H c , , 

PHQ = PH r Q; QHP = QH r P { ' 

QHQ = QH Q + QH r Q + QH X Q 

ultimately enters the answer in an even degree. So we arrive to the effective 
Hamiltonian acting in the model configuration subspace with the fixed num- 
ber of electrons in the d- and /-systems of a PTMC, but taking into account 
the LMCT and MLCT configurations in the «outer» configuration subspace. 

The resolvent eq. flBj) can be rewritten as a series with use of the bare 
resolvent: 

Ro(E) = (EQ-QHoQ)- 1 (8) 

derived from the operator QHqQ (conserving the quantities n<f and ni which 
are thus «good» quantum numbers). The first terms in the series yield: 

R(E) = (EQ - QH Q - QH r Q - QHxQ)' 1 « 
« Ro(E) + Rq(E) (QH r Q + QH X Q) R (E)+ (9) 
+R (E) (QH r Q + QH X Q) R (E) (QH r Q + QH X Q) R (E) 

When inserted in eq. (jSJ) the above expression corresponds to the fourth order 
of perturbation theory with respect to the resonance operator H r and at least 
second order with respect to H x . It is shown in Appendix [B] that the operator 
H r enters the effective Hamiltonian eq. (jSJ) in a twofold way. First, the H^ 
component of H rr which results in no ionization of the /-system (formal 
definition of the compoments HrT see Appendix [B]) acts within the model 
subspace ImP of our interest. Second, the H rr components of H rr take into 
account the «doubly ionic» MLCT and LMCT configurations in which the 
/-system is doubly ionized in either sense. We denote these configurations as 
(LM) 2 CT or (ML) 2 CT with respect to the direction of the charge transfer. 
These configurations have been shown to be important for estimating the 
amount of effective antiferromagnetic interaction. [42J Collecting the relevant 
terms yields: 
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H cS {E) « PH P + PH x P + PHg?P + 

+ Pi/Wg( 1+ )4 1+) g( 1+ ) J ff,g( 1+ )4 1+) g( 1+ ) J ff( + )p + 

+ PH^Q^Rt ) Q (1 - ) H x Q^Rt ) Q il - ) Ht ) P+, (10) 

+ PH^Q^R^Q^H^P + 

+ p^ r 2 -)g( 2 -)i?( 2 -)g( 2 -^( r 2 -)p 



which represents the approximate effective Hamiltonian acting in the sub- 
space with the fixed distrubution of electrons between the d- and /-systems. 
For derivation of its individual terms see Appendix iBl 

2.2.2 Magnetic limit 

In the previous Subsection we arrived to the effective Hamiltonian describing 
a PTMC in the configuration subspace with fixed distribution of electrons be- 
tween the noncorrelated (diamagnetic) ligands and the correlated (potentially 
magnetic) d-shells of the TMIs. The purpose of this and subsequent Subsec- 
tions is to derive the description of a PTMC in terms of its (/-shells only. This 
is analogous to the description of a mononuclear TMC by its <i-shell only, 
which is precisely the phenomenological crystal field theoryJTJ [2TJ [22] as se- 
quentially derived in the frames of the EHCF approach. [23j Like there the 
variables referring to the /-system must be integrated out by taking an aver- 
age over the approximate ground state \<&l) of the /-system. The only formal 
requirement for the wave function is that it is a single determinant one: 



It can be thus either closed shell RHF, ROHF, or the UHF wave functions. 
The subscripts A, a refer respectively for the /-MO and the spin projection. 

The uncertainty in defining \<&l) represents a problem. In the mononu- 
clear TMC's the definition of the configuration subspace in terms of the 
number of (/-electrons rid was sufficient to define the effective electrostatic 
field induced by the TMI in the ligands, so that the ligands' ground state 
could be calculated in this field with sufficient precision. In a PTMC the 
situation at this point is not that clear since electrons are allowed to travel 
between the d-shells of different TMI's and such a redistribution produces 




(11) 
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too strong perturbation of the /-system to be accomodated by a single single- 
determinant wave function. For that reason we have to make an additional 
assumption, conforming, however, to our purpose of describing magnetic in- 
teractions. We assume that there exists a preferred distribution {nd(i)} of 
electrons among the TMI's in the PTMC Fixing this distribution defines 
the Coulomb field felt by the electrons in the ligands and thus \<&l) can be 
calculated. 

To formalize this move, which in fact reduces the description of the PTMC 
to the model configuration subspace where only the local electronic motions 
are allowed we define the projection operator acting in the subspace of the 
n^-electronic states spanned by all possible determinants formed from AOs 
of the (i-system with the fixed distribution of c/-electrons: 

<$ = 0®Ii, (12) 

where // is the identity operator acting in the space of configurations of the 
/-system and O projects the configurations of the d-system on the subspace 
with the fixed numbers of electrons in each of the d-shells as formally defined 
by eq. (1421) in SectionO The operator 13 is complementary to : ^3+0 = l. 

With so defined projection operators we perform a further Lowdin parti- 
tion as applied to the effective Hamiltonian H eS (E) eq. ([5]) which results in 
another effective Hamiltonian: 

S)(E, <E) = <P# eff (E) ¥ + ¥H cS (E) Qm(<E)£lH cS (E) *p (13) 
The resolvent is naturally of the form: 

<K(<s) = [ae-OH^d]- 1 (i4) 

The action of the projection operators and ^3 on the terms of the effective 
Hamiltonian H eG (E) eq. (j3J) is: 

ypH py + ypH x py ^ 

yPH rr Py ^ 

ySPHoPQ. = QPH Py = (15) 

qSPHrrPQ. ^ 0; QPH rr Py ^ 

QH cS £l = 0P {Hq + H„ + H x ) P0 ^ 

Only the nonvanishing term ^PH rr PQ. and its hermitean conjugate admix 
the states with charge transfer between the d-shells i.e. metal to metal charge 
transfer states (hereinafter abbreviated as MMCT) laying in the «outer» 
configuration space ImPH to the states in our «model» space ImP*}3. 
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2.2.3 Ligand integration and ligand polarization 

Now we are almost ready to perform the final moves to eliminate the /- 
system variables. In variance with our original way[23j of reducing the ef- 
fective Hamiltonian of the type of eqs. (j3J) and (fT3|) acting in the restricted 
configuration space to one describing the mononuclear TMC in terms of its 
only d-shell which led us to the formal equivalent of the crystal field theory 
accompanied by the reliable estimates of the splitting parameters induced 
by this field|23j the situation in the PTMC's is more complicated since it 
requires more complete dressing procedure for the (/-electrons which even- 
tually includes the effects of polarization of the ligand sphere (diamagnetic 
media in which the TMI's are immerced) which accompany the virtual elec- 
tron transfers from one <i-shell to another. Such processes did not exist in 
the mononuclear case and for that reason had been originally neglected. It 
has been, however, explicitly shown that even in the mononuclear setting the 
polarization effects are responsible for the known nephelauxetic effect, (HI HE] 
which is the renormalization of the intrashell Coulomb (i.e. two-electronic) 
matrix elements due to indirect interaction between the fluctuations of the 
electronic density in the d-shell through the polarizations they induce in the 
ligands.[50j One may be sure that the partition described by eq. (jSJ) does not 
lead to the description of the TMC by the wave function of the form eq. ([3]) 
rather by a description through a superposition of the functions of that form 
each corresponding to fixed number of electrons in the subsystems, but still 
entangled. In fact the operator P projects out the states with the charge 
transfer between the d- and /-system, but does not assure the functional 
product form of eq. ([3]). Products of different states of the d- and /-systems 
satisfying only the condition of the fixed number of the electrons in each of 
the singled out groups enter in the expansion of the true ground state of the 
effective Hamiltonian eq. (J3J). The form eq. ([3]) is a kind of self consistent 
field approximation to it.[24j In order to improve this description one has to 
perform one more projection namely to the subspace of the products where 
all possible states of the /-system are replaced by its ground state. This is 
done as previously with use of the Lowdin projection using further projection 
operators. Let V be the projection operator on the ground state of the 
/-system: 

V = I d ® (16) 

complemented by the Kronecker product with 1$ - the identity operator with 
respect to the variables of the d-system. The complementary projection 
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operator is Q = Pty — V . Then applying the Lowdin partition as in the 
previous Subsections we arrive to the effective operator 1-L e ^{E, <£, Q) which 
acts in the subspace of the configurations where the /-system is fixed in its 
ground state 

H cS (E,£,n) = Vfi(£)V + VS)(£)QH(n)Qf)(£)V. (17) 

The Hamiltonian eq. (TTTj) does not affect the variables in the /-system as 
described by its ground state Thus the effective Hamiltonian for the d- 

system can be easily obtained by dropping the bra- ($>l\ on the right and the 
ket on the left in eq. (fTTl) (namely in this order). Then the remaining 
averages over the wave function yield the sought effective Hamiltonian 
for the d-system only (the variables related to the /-system are integrated 
out). In the subsequent Sections we identfy in turn the individual contribu- 
tions to it although the precise derivation of corresponding expressions will 
be kept in the respective Appendices. Further moves are stipulated by the 
projection to the single determinant ground state of the /-system as given by 
eq. f)17p . With use of the above generalized operators one can easily write 
the excited states relative to the single determinant ground state. 

For completeness we calculate the resolvent part of the operator eq. (fT7|) . 
We notice that for the ground state of the vanishing total spin (closed shell) 
assumed hereinafter the following holds: 



VH X V = 0; 
VH X Q + 

and calculate the second expression (or hermitean conjugate to it). It is also 
true that VHr? Q ^ and thus the operator Hrr contributes to the second 
order according to 

TZs(£l) and TZt{Q) contain respectively the projection operators to the singlet 
and triplet subspaces in the ImQ subspace. The average ( 



contributes to the effective Hamiltonian for the <i-electrons only. 
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3 Specific contributions to the c?-system effec- 
tive Hamiltonian. 



As shown in the previous Section the contributions to the effective Hamilto- 
nian for the G?-system of a PTMC stem from averaging eq. ( ITT)) over As 
we mentioned in the previous Subsection the leading terms are those which 
appear from the first term of eq. ( TTTT) or equivalently the average of the entire 
Hamiltonian eq. (TT3|) . In the subsequent Subsections we consider them one 
by one. 



3.1 Effective crystal field 

The simplest contribution comes from the averaging the first term in eq. ( TT3i) 
directly contributing to the first term of eq. (TTTj) over |$l).[23] The latter in 
its turn contains two contribution: the average of the Coulomb interaction 
over the ground state of the /-system \$>l) which trivially (see Appendix 
ID.1.3P yields the sum of the ionic/intraatomic contributions to the effective 
crystal field felt by the d-shells of the individual TMIs in a PTMC. It is 
consistently shown that the ionic contribution does not exceed 10% of the 
observed splitting. [151 HE ETJ EE E2] 

Nontrivial (and eventually the dominant) part of the latter is given by the 
average of the term H$ of the second order in H r originally performed while 
deriving the EHCF method for mononuclear TMCs.[23] Those derivation is 
generalized by treating individually the terms diagonal with respect to the 
TMI's indices ij as shown in Appendix ID. 1.41 

($ L | wh^pw \$ L ) = Ej V6f(j) = E,- E E * , «. 

{n K DW(jK)dt T (j)d, T (j) - (l-n K )D(^(jK)d+ T (j)d„ T (j)} 

This operator due to projection operator represents precisely as denoted 
the sum of the covalent contributions to the effecitve crystal fields felt by the 
d-shells of the individual TMIs in the PTMC. Then the overall result of the 
averaging of eff *P reads: 

($ L \yP(H + H rr + H x )P¥\$ L ) 
= Zi{H d (i) + V™(i)+VSl(i)+V^(i) (19) 

-2E^M<A(«)(*i|^(0|^)} 
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The summation in the last term extends to the d-AOs \i and /-AOs I of the 
i-th TMI. It appears only if the spin density on the /-th AO of the z-th atom 
is nonvanishing. It is possible for the ROHF and UHF functions \<&l) only 
and it is new as compared to the original EHCF[23] where only the RHF form 
of was allowed. However, hereinafter we concentrate on the systems in 
which in the zero approximation the ligands are diamagnetic, i.e. have no 
magnetic moments/spin densities in the /-system (and are thus representable 
by the RHF function) so that this term further is not going to reappear. 

3.2 Magnetic contributions 

In the previous Subsection we rederived the second order contribution of the 
one-electron hopping between the d- and /-systems to the effective Hamil- 
tonian. Its part which is diagonal with respect to the TMI's indices % and 
j (or, equivalently, one acting within the configuration subspace Irn^P) re- 
sulted as expected in the sum of the crystal field contributions for individual 
TMIs. The r/'-off-diagonal part of Hr? (that according to the formal expres- 
sion tyHrr^O. coupling the configurations in the subspaces Im^3 and ImH) 
contributes to the effective operator for the d-system according to: 



by which the MMCT configurations are projected out. Formally it is nothing 
but the average of the second term in eq. (fT3l) . Eventually it represents 
the leading magnetic contribution. Details of the derivation are given in 
Appendix ID. 1L The corresponding configurations are shown in Figs. [TJ - HJ 
They describe the situation which is conveniently described as effective one- 
electron transfers between different d-shells. However, it has been shown[42j 
that other processes having nothing to do with transfers of electrons between 
interacting (i-shells rather those which can be described as correlated electron 
transfers between the /-system and the d-shells in which two electrons of 
the opposite spin projection are excited from a single /-MO to two (i-shells 
under consideration. The correlations in this case is as well of «kinematic» 
nature and thus they contribute predominantly antiferromagnetically. The 
corresponding configurations are shown in Fig. [5j these contributions appear 
as a result of averaging the terms of the second order with respect to the 
ionizing components Hr^ 1 in eq. (|40p of Appendix [Bl 



( 




) 



(20) 



15 



The terms eqs. (1561) . (I58p must be summed with the spin dependent fac- 
tors assembled in Table 0bf Appendix ID. 21 and accroding to one of the four 
"cases" (i) - (iv) identified by Weihe and Gudel|4Ti] to take into account the 
effect of the different variations of the spins in the (i-shells involved in the 
electron transfer processes up on the contributions to the effective exchange 
between these d-shells. As one can see from Table H] the terms with different 
variation of the local spins contribute with different sign. The classification 
to the "cases" (i) - (iv) is based on the assumption that the transfer of an 
electron involving a half-filled orbital always result in a reduction of the local 
spins by 1/2 in either of the involved ^-shells irrespective to the sense of the 
transfer. If either of the involved orbitals is either completely occupied or 
empty the electron transfer process may result either in increase or decrease 
of the local spin by 1/2 in the d-shell where such an orbital occurs. This 
results in only one contribution for the case (i), two contributions for the 
cases (ii) and (iii), and four contributions for the case (iv). These contribu- 
tions enter with numerical unpairity factors specific for each case and with 
the signs specific for each combination of the possible variations of the local 
spins. The values of the energy denominators (D^ liz \ D^ 2± \ T)^) referring to 
the intermediate states in the outer configuration subspaces must be taken 
according to the variation of the local spins specific for each contribution 
(that is to be higher by the intrashell exchange energy when the spin goes 
down and to be smaller by the same quantity when the spin goes up). That 
structure of the contributions and of their combinations entering the final 
expressions for the exchange parameters moved the authors[42j to expand 
the denominators' products against the intraatomic exchange energies. Due 
to alternating signs of the contributions corresponding to the different conbi- 
nations of the intermediate local spins the terms of the lowest nonvanishing 
order survive for each of the cases (i) - (iv) which are respectively of the 
zeroth power for the case (i), of power one for the cases (ii) and (iii) and of 
the second power for the case (iv). This is how the expressions for the cases 
(i) - (iv) have been derived previously |42| We do not use these expansions 
and employ "exact" values of the energy denominators. 

Our derivation allows an additional classification according to the types 
of the terms in the sum over k\ which run through the /-MOs. These types 
(a) - (d) formally defined by eq. (157|) correspond to the four possible transfer 
paths shown in Figs. [1] - H] and respectively involving two occupied, two 
empty, and one empty and one occupied /-MOs in various orders. They 
describe contributions of the states in the MMCT configutation subspace 
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Table 1: Generalized factors of the matrix elements contributing to the ex- 
change parameters. 



a 


n K n x D^ + \j\)D^ + \jK) 


b 


-n K (l-n x )D^(jK)D^(z\) 


c 


- (1 - n K ) n x D^-\iK)D^ l+ \]X) 


d 


(1 - nx) (1 - n K ) DV-*(iK)DV-)(i\) 


e 


n K n x DW (m) + (jX) 


f 


(1 - n x ) (1 - n K ) D (1 -\i\) DV~\iK) + D^(jX) 



and thus contain the energy denominator S^,-. 

The values of these factors are generalized and summarized in Table [TJ 
Each of these matrix elements contains the following expression: 

Mi) /WO /WO A/A 0?) 

as a common factor. The first term in square brackets corresponds to an 
electron transfer from the z'-th d-shell to the j-th one. The second term 
comes from the process of spin-correlated transfer of two electrons with op- 
posite spin projections from the occupied Z-MOs to the two d-shells by this 
effectuating coupling between them. 

In the previous Sections we were able to basically rederive and generalize 
the perturbative formulae for the contributions to the exchange in a form 
suitable for programming. They are also complemented by additional terms 
stemming from various processes involving polarization of the /-system. As 
expected the terms are numerous. However, additional considerations allow 
to range these terms in order of their importance. At the present first step 
we concentrate on the estimates of the contributions stemming from the 
electron transfers which are formally of the fourth power with respect to the 
one-electron hopping integrals /3's. 

4 Details of implementation & calculation re- 
sults 

The formulae for the effective exchange parameters derived in the previous 
Sections do not require any additional quantities except those which are 
already calculated in the context of the EHCF method. |23j These are the 



^i-tj ^ "'k'I'X-^ ' K x->ij 



(2+) 



(21) 
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Table 2: Calculated and experimental values of the exchange parameters in 
/x-oxo bridged Cr(III) dimers. 



Compound (CCSD Code) 


J exp , cm 1 Ref. |60J 


Jcaici cm 


Geometry source 


[(NH 3 )5CrOCr(NH 3 ) 5 ] 4+ 


450 


408 


[52J 


GAMTUJ 


510 


424 


[53J 


VIDTIL 


100 


190 


[54J 


ZUVMIM 


60 


150 


[55J 



/-MO LCAO expansion coefficients, respective orbital energies, d-AO-l-MO 
resonance integrals Pfj,\(i), etc. We implemented the derived formulae as a 
program suite accepting standard quantum chemical input (molecular com- 
position and geometry) using the GEOMO package[51J (QCPE No 290) as 
a source of the subroutines for performing the calulations of the molecu- 
lar integrals and performing semi-empirical SCF MO LCAO procedures for 
the /-system. The package has been tested agaist the compounds of the 
Cr(III)OCr(III) family[52j [53j [541 [55] known as «basic rhodo» compounds 
which have been synthetized by S.M. Jorgensen 130 years ago.[56j The pa- 
rameterization procedure follows general EHCF methodology and will be 
described elsewhere. 

We performed a series of calculations for those compounds studied in 
the phenomenological setting where the structural data were available. The 
results relative to the /^-oxo bridged Cr(III) dimers are given in Table [2J 

As one can see the order of magnitude of the exchange parameters in 
this series of compounds is correctly reproduced as are their trends de- 
pendent on the chemical composition and bridge geometry (basic erhythro- 
compound [(NH3)5CrOCr(NHs)5] 4+ and GAMTUJ have linear bridge geom- 
etry whereas VIDTIL and ZUVMIM are bent with the CrOCr angle in the 
range 128°-^132°) , although the amplitude of the angular dependence as 
obtained in the calculation is somewhat smaller than the experiemental one. 
Analogous calculations performed for the /z-oxo bridged Fe(III) dimers (Ta- 
ble [3]) show similarly reasonable agreement between the experimental and 
calculated values of the exchange constants. 

The most remarkable feature we addressed in this context is the effect of 
protonation of the oxo-bridge upon the magnitude of the exchange constant. 
Comparing the calculated values in the last two lines of Table [3] we see that 
the protonation as expected breaks at least one of the superexchange paths 
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Table 3: Calculated and experimental values of the exchange parameters in 
/x-oxo bridged Fe(III) dimers. 



Compound (CCSD Code) 


J cxp , cm" 1 Ref. [41J 


Jcaici cm 


Geometry source 


DIBXAN 


238 


240 


[62J 


VABMUG 


264 


256 


164] 


PYCXFE 


214 


118 


[61J 


COCJIN deprotonated 


242 


110 


[63J 


COCJIN 




30 


[65J 



going through the oxo-bridge which as expected as well significantly reduces 
the magnitude of the effective exchange parameter. 



5 Conclusion 

By this we extend our work [23J of the year 1992 devoted to calculating 
the intrashell excitations in the d-shells of coordination compounds of the 
first transition metal row, which resulted in the Effective Hamiltonian Crys- 
tal Field (EHCF) method, to their polinuclear analogs in order to assure 
the description of several open d-shells and of magnetic interactions of the 
effective spins residing in these shells. This is a challenging task since it 
requires improving the precision of ca. 1000 cm -1 (that of describing the 
excitation energies of the single rf-shells by the already well successful EHCF 
method) to the that of ca. 10 -j- 100 cm -1 characteristic for the energies 
required to reorient the spins i.e. eventually by two orders of magnitude. 
This is performed within the same paradigm as used for the EHCF method: 
the concerted usage of the McWeeny's group-function approximation and 
the Lowdin partition technique. These are used to develop the effective de- 
scription of the (/-system of the polynuclear complexes composed of several 
d-shells including the working formulae for the exchange parameters between 
the (/-shells belonging to different transition metal ions. These formulae are 
implemented in the package MagAixTic and tested against a series of binu- 
clear complexes of trivalent cations featuring /i-oxygen superexchange paths 
in order to confirm the reproducibility of the trends in the series of values 
of exchange parameters as well as the magnitude of these parameters for 
the compounds differing by the chemical substitutions and other details of 
composition and structure. The results of calculations are in a reasonable 
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agreement with available experimental data and other theoretical methods. 
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A Hamiltonian contributions 



A.l Subsystem's Hamiltonians 
A. 1.1 Hamiltonian for the d-System 

The bare Hamiltonian for the (i-system reads: 

E^) + ^E y coui(y) 

The Hamiltonians for the individual (f-shells are taken in the atomic sym- 
metric approximation dating back to Refs. |66] W0\ |69| |68| 167] so that the 
Coulomb interaction and the exchange splitting of d-electrons are decribed 
by two parameters (U dd (i) is the one-electron core-attraction parameter): 

H d {i) = U dd (i)n d (i) + ^g dd (i)h d (i) (n d (i) - 1) - K M {i)S 2 d (i) (22) 

where h d (i) is the operator of the number of particles, S d (i) is the operator 
of the total spin, U dd (i) is the one-electron core attraction parameter, g dd (i) 
is the average parameter of the Coulomb interaction, K dd (i) is the average 
parameter of electronic exchange, all referring to the i-th (f-shell. The main 
advantage of this form of the Hamiltonian is that it preserves a higher SO (5) 
rather than the actual SO(3) symmetry of the atom so that it remains in- 
variant under arbitrary 5x5 orthogonal transformation of the one-electron 
d-states. This will be necessary when going to the basis of the eigenstates of 
the local crystal field operator. 

With such a bare Hamiltonian the state of the i-th d-shell is uniquely 
characterized by two quantum numbers: n d {i) - number of electrons, we use 
to build the configuration subspaces and to decompose the system in parts, 
and S d (i) - the total spin of the given <i-shell. Then the energy of the entire 
manifold of the states with given n d (i) and S d (i) is given by: 

U M (i)n d (i) + ^g M (i)n d (i) (n d (i) - 1) - K M (i)S d (i) (S d (t) + 1) (23) 

This form assures for the Hund's rule for the states of the (f-shells due to 
positiveness of the intrashell exchange parameters K dd (i) and for the com- 
pliance of the bare spectrum with the Lande interval's rule as suggested in 
Refs. |15| W7\ . The values of the core attraction parameters are taken as 
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implemented in the EHCF package |23j . Those of the average interaction 
parameters are given by J0rgensen as well: 

9dd {i) = F° dd (i) - 1 (F 2 M (i) + *&(«)) 

K dd (i) = § B(i) + ~C(z) = A (**(<) + F^{i)) 

- not a very great improvement over Refs. [39, 40, 431H2] where B(i) and C(i) 
are the Racah parameters for the z'-th TMI and Ff v (i) are the corresponding 
Slater-Condon parameters. 

A. 1.2 Hamiltonian for the /-system 

INDO parameterization for the first row elements has been introduced in 
Ref. [72] . Its extensions to the transition metal atoms had been proposed 
in Refs. [73| IT4"| [75] . Within that setting the main problem was to im- 
plement intraatomic Coulomb and exchange two-electron integrals (Slater- 
Condon parameters F k , G h ) since wihtin the INDO (SCF) setting they cover 
all necessary intraatomic parameters allowed by symmetry. Rinaldi[51J have 
shown that the intraatomic configuration interaction involving s-, p-, and d- 
(sub)shells requires additional intraatomic integrals. These latter, however, 
again disappear if the ^-shells enter the wave function as direct multipliers 
eq. ([2]). This is thus the case for the present method. 

Then in the second quantization form the Hamiltonian for the /-system 
reads: _^ 

m,<7 L 

+ E E {Uu - £ v w -v LM )iX+ 

L leL,cr L'^L (24) 

+ E 0mi(m+l a + h.c.) + £ MX+ 

ml, a IV ', a 

+5E«v(" / |i l, nW«; 

here (l a ) are the creation (annihilation) of an electron with the spin pro- 
jection cr on an /-AO. First term in the expression eq. (j24"|) describes the 
interaction of the 4s- and 4p-electrons of the metal (m = 4s, 4p x , 4p y , 4p z ) 
with the metal core (parameters U mm < 0) and the ligand atoms cores (pa- 
rameters Vml > 0). Second term describes interaction of the ligand electrons 
with the ligand cores (parameters Uu < 0), with the cores of the other ligand 
atoms (parameters Vlu > 0) and with the metal core (parameter Vlm > 0). 
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Third and fourth terms describe the resonance interactions in the ligand sub- 
system (parameters f3 m i < and (3w < 0). Last term describes the Coulomb 
interactions between electrons ((//' | l"V") are the corresponding two-electron 
integrals) . 

A. 1.3 Fockian for the /-system 

The calculation of the wave function of the /-system is the prerequisite 

Variation principle applied to the effective Hamiltonian with the trial 
function of the form eq. leads to the self-consistent system of equations: 

H e / f ® d = E d § d 

a n n n 

H?"& = E l ¥ 
H e / f = H d + | H c + H rr | $(,) , 1 ; 

In the above system the effective Hamiltonian ti e J for the rf-electron sub- 
system depends on the wave function of the ligand subsystem $q, and in 
its turn the effective Hamiltonian i/j for the ligand susbsytem depends 
on the (i-electrons' wave functions $q. These equations must be solved self- 
consistently as well. In the EHCF method Ref. [22] the Slater determinant 
$o, is constructed of MO's of the /-system, obtained from the Hartree-Fock 
equations in the INDO approximation for the valence electrons of the ligands. 
In this case the transition from the bare Hamiltonian for the /-system to 
the corresponding effective (dressed) Hamiltonian reduces to renormalization 
of one-electron parameters related to the TMI: 

eff " (26) 

where Ua is the parameter of the interaction of 4s and 4p-electrons (i = 4s, 
4p x , 4p y , 4p z ) with the TMI core, Z M is TMI core charge, = (fifi\ii) — 
are the parameters of intraatomic Coulomb interactions. The $q 
function thus obtained is used further for constructing the effective Hamil- 
tonian for the (i-shell. 

A. 2 Intersubsystem interaction operators 

Here we introduce explicit definitions for the interaction operators acting 
between the d- and /-system of a PTMC. 
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A. 2.1 dZ-Resonance operator 

The resonance operator H r describing one-electron hopping between the d- 
shells and the ligands has the form: 

H r = + 
H r +) = -EEE^WCWa*,, (27) 

a fi,i X 
c fx,i X 

where A runs over /-MOs, P^i) is the resonance (hopping) integral between 
the /i-th d-AO of the i-th TMI and the A-th /-MO, the fermion operator 
d+(i) creates an electron with the spin projection a on the /i-th d-AO of the 
i-th TMI, and a\ a creates an electron with the spin projection a on the A-th 
/-MO. One can easily check that the terms are the hermitean conjugates 
of each other. Action of apparently results in positive ionization of the 
/-system (/ — >■ d one-electron transfer), that of Hr ^ refers to the negative 
ionization of the /-system (d — > I one-electron transfer) . The above definition 
can be somewhat simplified by using the spinor notation: 

H r +) = -EE'WWOfcA, ( 28 ) 

fi,i X 

H r~ ] = -EE'WOtfM/.w- 

H,i X 

A. 2. 2 d/-Coulomb and exchange operators 

Although fundamentally we rely upon the INDO approximation for the bare 
Hamiltonian for the PTMC we have to make certain concessions and regroup- 
ing of terms in order to profit from the symmetries characteristic for the 
atomic problem. Specifically we separate the Coulomb interaction into sym- 
metric (superscript "(s)") and asymmetric (superscript "(a)") parts of which 
the first incures only the uniform shift of the d-\eve\s in each given TMI 
whereas the second induces the splitting of the otherwise degenerate cZ-levels. 
The symmetry mean by the symmetric part is that of the SO (5) - i.e. of an 
arbitarary orthogonal transformation of the d-orbitals. These two contri- 
butions further subdivide into an interaatomic (marked by subscript «1») 
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part describing the interactions with the electrons in the noncorrelated sp- 
AOs of a given TMI and the interatomic one describing the interactions with 
electrons residing on other atoms (marked by subscript «2»). With these 
assumptions the Coulomb interaction between the subsystems of a PTMC 
reads: 



i L 

H£ = e^^n^rUFo^zi) (29) 

i L 
ifi L 



mm' fi 

where n^ii) and hi are the operators of the numbers of electrons in the i-th 
(i-shell and on the Z-th /-AO; h^{i) is that for the /z-th d-AO of the i-th TMI; 
L stands for all orbitals I centered on the ligand atom L; 7 = x,y,z g^i) are 
the energy parameters (in fact - the average Coulomb interaction integrals) 
and the interatomic contributions V mm r are given by: 

Knm' = ^ ] ^ ] ^mm'^mcr^rn'^r W 

with 

= F>\/&TT \ , 

VW+i V 



L lei 

a.r 



V£,=*/S liDH J_ m L )F i iR L )Yr"'(0 LM . 



k=2A 



(30) 

where (Rl,&l,4>l) are the spherical coordinates of the ligand atom L (rela- 
tive to a TMI consequently located in the center of the coordinate frame); 
ym-m (# L; L ) are the spherical functions with the phases defined following 
Condon and Shortley. [7T| Functions F^Rl) are the integrals of squares of 
the radial parts R n i(r) of the atomic d- functions: 

R 00 

F k (R) = iT (fc+1) j r k R 2 nl (r)r 2 dr + R k j r^ k+1) R 2 nl (r)r 2 dr (31) 

R 
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depend on the distance from the atom of metal to the atom L. For the 
Slater AOs R n i(r) the functions Fk{R) are explicitly known. |22j 

Experience of usage of the EHCF numerical procedure acquired so far 
shows that the effect of the Coulomb splitting (asymmetric part of the 
Coulomb interaction) can be safely neglected when it goes about estimat- 
ing the energy denominators in the resolvent in Section IC.ll which is done 
consistently. 

The exchange operator decomposes analogously with that difference that 
in the INDO approximation only the intraatomic contribution on the TMIs 
is present. 

H x = Hjp+HW (32) 

We however restrict ourselves by the symmetric part of the intraatomic in- 
tershell exchange for the time being: 



ifW = -2J2 Yl \K sd (i)S s (i) + K pd (i)S p (i)] (33) 

i fj,ed-AO 
iGl-AO 

s d (i) = x;^(o 



where S^i) and Si(i) are the electron spin operators for electrons in the /i- 
th d-AO of the i-th TMI and the Z-th /-AO of the i-th TMI. The exchange 
parameters are given by: 

K sd {i) = \G 2 sd {i) 

9 113 

K pd (i) = -Gl d (i) + Y^Gl d (i) (34) 

B Expansion of the resolvent eq. (Ej) 

Inserting the expansion of the resolvent eq. in eq. (j3J) we obtain: 

H cS (E) = PH P + PH X P + 
+PH r QR QH r P + 

+PH r QR QH x QR QH r P + (35) 
-\-PH r QRoQH r QRoQH r QRoQH r P + 
+PH r QR QH x QR QH x QR QH r P + ... 
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In the above expression (and in its continuation whatsoever) the odd powers 
of H r disappear since the operator H r changes the number of electrons in 
the d- and /-systems by one in either sum, and in case when the number of 
such changes is odd the total number of electrons in the d- and/or /-systems 
cannot be conserved. Thus the number of the H r multipliers must be even. 

Further analysis can be based on the notion of "graduality" of the comple- 
mentary projection operator Q and the resolvent expressed by the following 
expansion for these quantities: 

Q = 0Q (n±) 

Ro = Q R (n±) ( 36 ) 

n 

where n refers to the number of electrons added to ("+") or taken from ("— ") 
the cZ-system (this notation precisely refers to the degree of ionization of the 
/-system), and © stands for the direct (block) sum of the corresponding 
matrices. For the terms in eq. (136|) the following orthogonality conditions 
hold: 

Q (n±) R (m±) Q (p±) = 6nm6mpR (m±) 
Q(m±) R (n±)p = p R (n±)Q(m±) = q 

With this decomposition for the resolvent and graduation of the ImQ 
configuration subspace the effective Hamiltonian rewrites: 

H cS (E) = P(H + H X )P + 

+ PH r Q^R^Q^H r P + 

+ PH r Q {1±) R (1±) Q (1±) H X Q {1±) R (1±) Q {1±) H r P + 

+ PF r (5 (1±) .R (1±) Q (1±) H r Q^ 2±) i? (2±) Q (2±) F r (5 (1±) i? (1±) Q (1±) H r lftf) 

+ PH r QR QH x QR QH x QRoQH r P + ... 

where the terms with «+» and «-» superscripts are summed separately. 

As for the terms containing the powers of the exchange operator H x QRq 
only they seem to be summable since they involve only the terms of the same 
graduality. It looks like: 
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PH r QR QH x QR QH r P + 
PH r QRQQH x QR QH x QRQQH r P + 
PH r QRoQH x QRoQH x QRoQH x QRoQH r P + ... 



oo 



PKrQRo (QH x QR ) m QH r P 




which represents nothing but multiple scattering of an electron wandering in 
the 1-system by all possible TMIs in all possible combinations. For the time 
being we retain only the first order term in H x . Thus we obtain: 



+ PH r Q^R^Q^H r Q^R^Q^H r Q^R^Q^H r P + 

+ PH r Q^R^Q^H r Q^R^Q^H r Q^R^Q^H r p 



Since the model configuration subspace is that of the nonionized /-system 
the resonance operator eq. (j27p can enter in the answer in even powers. 
We classify the terms in the effective interaction stemming from the parti- 
tion procedure eq. according to the degree of the ionization these terms 
introduce to the /-system. In the lowest order we can write: 



H ea (E) 



« PH P + PH X P + 
+ PH r QR QH r P + 
+ PH r QR QH x QR QH r P + 



(38) 



H rr = H r R H r = + + 



(39) 



The ionization components of the operator H rr have the form: 



H?r +) = E E 0,x(i)^ K (j)d+(i)a Xa RWdt T (j)a, 



+ 



#w= e E M*)/Wi)x 



(40) 



{a+ d^(z)i2< 1+ >d+ UWr + 
+ d+ a {i)a Xa R^at T d VT {j)} 
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Comparing eq. (l3"Tj) with the definitions eqs. (1391) . (1401) we can identify 
the products P r g (1±) P (1±) g (1±) P r with components h£ ±] so that eq. (jggjl 
rewrites: 



H eS (E) « PH P + PH X P + PH^P + 

+ ph^q^r^ +) q^h x q^r^ +) q^hMp + 

+ p^-)g( 1 -)4 1 - ) g( 1 -)^g( 1 -)4 1 - ) g( 1 -)iyHp+ (4i) 

+ pp( r 2+ )g( 2+ )p( 2+ )g( 2+ )p( T 2+ )p + 

+ pp( 2 -)g( 2 -)p( 2 -)g( 2 -)p( r 2 -)p 

This represents the approximate effective Hamiltonian acting in the subspace 
with the fixed distribution of electrons between the d- and /-systems. 



C Resolvents and energy estimates in the outer 
subspaces 

C.l Resolvents for the configuration space partition 

The components for the resolvent acting in the configuration subspace Img 
for the subspace with the singly ionized ligands have the form: 

(we remind that the superscripts for the resolvents refer to the ionization 
degree of the /-system). They involve the projection operators in the d- 
system 

O = 0({n d (j)}); 
Q* = 0({n d (j)}^ k ,n d (k) + l); (42) 
O fc = 0({n d (k)}^ k ,n d (k)-l) 

are those projecting to the subspaces with the respective fixed numbers of 
electrons in the cZ-shell of the j-th TMI in the polynuclear complex. The 
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projection operators further subdivide into components corresponding to the 
resulting total spin of the affected d-shell, that is: 



t (s(k) + i) + o* (s{k) - . 



Due to cumbersome notation we shall not always indicate this subdivision, 
but shall keep track of it since it is necessary (see Sections |D] and [3]) • 

The projection and Fermi operators obey the following commutation 
relations: [23J 

d(j)& + = Od(j); d+(j)&_ = Od+U) u ^ 
& + d + (j) = d + (j)D; oLd(j) = d(j)Q [ ] 

The energy denominators: 

DV+)(kK) = [E - (I K - A{k) - g K (k))) 1 
D^(kK) = [E- (I(k) — A K — g K (k))\ - 1 

where g K (k) is the Coulomb interaction between an electron and a hole in 
the k-th (i-shell and K-th /-MO and the operators and J's and ^4's are the 
ionization potentals (IPs) and electronic affinities (EAs) of the ligands and 
d-shells. Here we as well will be keeping track of the spins of the k-th rf-shell 
resulting from the electron addition or abstraction. 

The estimates for the IPs and EAs for noncorrelated /-system MOs are 
taken according to the Koopmans theorem: 

I K ,A K = -e K (44) 

coming from the semiempirical SCF procedure as applied to the /-system. 
The ionization potentials and electron affinities for the d-shells (I(k),A(k)) 
are more complicated since they have to account for the difference between 
the electron addition and electron subtraction to/from the (i-shell as well as 
for the total spin of the (/-shell emerging as a result of either of the processes. 
Using the expression for the bare energies eq. (T2S|) for the relevant states we 
obtain: 

r /.N tt / n nu , U \ -n , / -K dd (k) (S d (k) + §) increasing spin 

I(k) = -U M {k)-g dd {k){n d {k)-l) + l /,w / M . n , 

[K dd (k) [S d (k) + jj decreasing spin 

A{k) = -^W-^(^ W + {^ (A;) ^ (A;) + ^ X inCTeaS ingSpin 

y-K dd (k) (S d (k) + decreasing spin 
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The difference between the IP's and EA's for the different resulting spins 
yields in both cases precisely ±JQd(/c) (2Sd{k) + 1) which with attention to 
the notation coincides with the result given in Ref. [ID]. By this the resolvent 
relevant for the partition eq. (j3J) is defined. 



C.2 Resolvents for the magnetic limit 



The resolvent appearing in the partition eq. (lTB"j) is the most peculiar one 
since as one can see from eq. (|47j) of Section [D] it describes the contribution 
of a process which can be called a dressed effective hopping and the relevant 
excited states in the configuration subspace ImH involve the states with dis- 
tribution of d-electrons different from the fixed one as well as the excitations 
of the /-system. The operators projecting on the basis states in ImH are 
thus: 



C— >A/ \ K— >X I 



K->A 



which are to be supplied with the respective energy denominators: 



(0) 

ijSK—t-X 



1 - 5 KX ) ( ~ Ar J) - 9xn + 2K XK 



(o) 

ijTK-*\ 



e - (/(/)- Mj)- gij + E^ 

1 - <U ( It 3) ~ At* ~ 9x 



(45) 



where and A(j) are the IP and EA for the i-th and j-th <i-shells respec- 
tively as defined in Section IC.ll They contain the dependence of the energy 
of the MMCT states on the spins of the d-shells involved in the transfer. The 
quantity is the average Coulomb interaction parameter for electron and 
hole residing in the respective d-shells, which can be taken in a simple form; 
g\ K and K\ K are, respectively, the Coulomb and exchange integrals for the 
pair of /-MOs coming from the semiempirical SCF calculation; is the 

reorganization energy of the /-system acquired under the transition, which 
can be estimated as: 



E R^ l] = J2^9iL~ 9jl) 
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where the effective charges qi in the /-system come from the semiempirical 
SCF calculation for the /-system in the electrostatic field induced by the fixed 
distribution of d-electrons. The ionization potentials and electron affinities 
Ik~* and are evaluated according to the Koopmans theorem as the 

orbital energies of the corresponding MOs shifted by the first order correction 
from the electrostatic field induced by the charge transferred between the d- 
shells. That is: 

L leL 

Further energy denominators required for taking into account the states in 
the configuration subspace ImQ( 2± )0 are those corresponding to double ion- 
ization of the /- and (/-systems: 

E - + -A K -A X + 9ij + g KX + 4^ J) ) 

1 - 5 kX ) (/<*-*> - - g XK + 2K XK )] 1 

E-(l K + I x - A{i) - A(j) + 9ij + g KX + 

1 - 5 KX ) (if^ - A[ 1 ^ - g XK + 2K XK )y 1 

The projection operators to the configurations in the doubly ionized subspace 
are: 



£,(2-) 



D 



(2+) 
ijSnX 



—>x 

T 7 (2- 

K— >X 



S(2+) 
>A 



T 7 (2-) 
>A 



C.3 Resolvents for the /-system polarization 

For the /-system only singly excited states over its singlet (closed shell) 
ground state are previewed. El [9] These states apparently classify accord- 
ing to their total spin which can be thus singlet and triplet. For that reason 
the resolvent eq. (1171) decomposes 



\®S \p(0) 



/jXk 

7 



V, 



(o) 

Tk-^X 



K->X 



(46) 
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where the superscript on the left side refers to the fact that the /-system is not 
ionized in this configuration subspace; |3>f_^) and 5>^w stand respectively 



for the singlet and triplet excited states of the /-system with one electron 
excited from the MO \k) occupied in the ground state to the MO |A) empty 
in the ground state and 7 denotes the component (x, y, z) of the triplet state. 
The excited states of the /-system yield the energy denominators 

VsLx = [n-{I K -A x -g XK + 2K XK )}- 1 
4°La = [n-ilK-Ax-g*)]- 1 

naturally independent on the projections in the triplet state. Here I K and A\ 
are respectively ionization potential and electron affinity of the /-system as 
defined by eq. (144)) . g\ R is the Coulomb integral for the pair of Z-MOs, K\ K 
is the exchange integral for the pair of Z-MOs. 

D Contributions to the effective Hamiltonian 

The contributions to the effective Hamiltonian for the <i-system as described 
in Section [3] are the averages over the ground state of the Z-system 
The averages of the diagonal (symmetric) parts of the interaction operators 
comes out trivially. Nontrivial contributions come from the operators capa- 
ble to excite the Z-system. We consider them one by one in the following 
Subsections. 

D.l Effective one-electron hopping between TMIs and 
related terms 

D.l.l Preliminaries 

We start with the action of the operator Hrr on the ground state of the 
Z-system. It yields: 

fj,iXcr ujKT 

wMO £ ®\D^XkK')a KlT ,\<b L ){<b L \at, T 4t T {3Wr\<$>L) 

kK'r' 



x < a 



kK'r' 
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Employing the anticommutation relations for the Fermi operators: 



and the values of the averages of the products of the latter over the single- 
determinant functions: 

a+, T ,a KT |$ L ) = S KK ,S TT/ n KT ; ($ L | a K , Tl a+ T \® L ) = 5 KK ,5 TT >{1 - n KT ) 

we get: 

fiiXa vjnr 

EMaA^Xt \$ L ) D^ + \kn)O k + d+ T (j)n KT + 

k 

~ E S kj d+ r (i)a Jur a+ r |$ L ) D^-\kK)O k _ d UT (j)(l - n K1 

k 

Complementing this by the projection operator O from the left and per- 
forming its commutation with the Fermi operators according to eq. (j43|) we 

get 

fiiXaujKT 

+ 0* + Oid+ (i)d VT {j)a Xa at T D^-){jK){\ - n K )} 

which describes the process which admixes the states in the configuration 
subspace luiPVpP to those in the subspace YmVQQ^ 1 ^. It is a starting 
point for further evaluations. 

Due to the definition of the projection operators *P and the Fermi 
operators referring to the /-systems commute with them so that they act 
directly on the wave function an d produce the singlet |$f_> A ) and three 

components of the triplet excited states of the ligands. This shows 

that in the general case a motion of an electron between the d-shells is 
accompanied by the charge and/or spin polarization of the ligands. It had 
been checked that the excitation energies in the /-system of the exemplary 
molecules considered in the present paper (see Section H|) can be estimated 
of about 10 eV, thus their contribution is strongly damped. 

An important exception is represented by the terms in eq. (|47p where 
K — A. In this case a = t also holds and the /-system incidentally acquires 
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no excitation. This singles out the following terms from eq. (T4T|) : 



E E/Wi)A*0')x 



i^j/xi/ At 



(48) 



{OiQP + d„ T (i)d+ T (j) ,D (1+) (jA)n A + 



This expression generalizes and provides the numerical estimate to the matrix 
elements responsible for the MMCT admixture. [33, SHI S2J EE] 

This is apparently the one-electron operator acting in the (i-system and 
describing electron transfers between the <i-shells of different TMIs through 
the /-system. Two groups of channels for such transfers do exist: the occupied 
and empty MOs of the /-system. In the first case (occupied n-th /-MO) an 
electron is captured from it by the j-th c/-shell and then the hole thus emerged 
is healed on account of an electron coming from the 2-th <i-shell. In the second 
case (empty /t-th /-MO) an electron leaves the j-th cZ-shell and then lands in 
the i-th c?-shell. Despite an unsymmetric appearance this effective operator 
is hermitean. The terms with % = j were precisely the covalent contributions 
to the effective crystal field. |23j 




D.1.2 Averaging over the /-system 

As in the original EHCF derivation[23j the transition from the effective 
Hamiltonian for the (P)TMC acting in the configuration subspace with the 
fixed distribution of electrons between the d- and /-systems to that for the 
(i-system only is performed by averaging the former over the ground state of 
the /-system represented by the single determinant wave function In 
variance with the case of TMC in PTMC we additionally assumed that the 
must be calculated for the fixed distribution of the ^-electrons among 
several d-shells of a PTMC. Formally it is reflected in the requirement that 
the effective Hamiltonian for the PTMC acts in the model configuration sub- 
space ImP*p. Introducing a shorthand notation 




«...», = (*l\.~\*l) 



we easily arrive to 
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D.1.3 Ionic contribution to the effective crystal field 

It is necessary to average the operators H c and H rr with the function of the 
ground state of the /-system <3>£. For ((H c ))i we get an expression: 

= EE g^Pu^ + E E v^p LL d+j ua . (49) 

The first term in this expression describes the shifts of the <i-levels coming 
from the interactions of c/-electrons with electrons on the 4s- and 4p-orbitals 
of the metal atom. The second term represents the interaction of <i-electrons 
with electrons on the valence orbitals of the ligands. The sum of the first and 
the second terms in the expression eq. () has the form of the operator of the 
crystal field induced by the effective charges located on the ligand atoms. 

D.1.4 Covalent contribution to the effective crystal field 

Multiplying eq. (|4T|) by ($£ | from the left and integrating over the variables 
pertaining to the /-system and taking into account the above averages over 
the Slater determinant wave functions: 5 K \5 TU n K and 5 K \5 T(T (1 — n K ) in the 
first and the second rows of eq. (147|) . respectively, we get (the occupation 
numbers are always unity or zero): 



($ L |/4 0) I$ L >= EEWM X 

fj,iuj KT 

{-OtOP+n.D^ij^d^d+^j) + ( 5 °) 
■ OVOL(l - n K )D(^(j K )d+ T (z)Mi)} 



for the terms with % ^ j and the following: 



fll/J KT 



{n K D^(jK)dt T (j)d,T(j) - (l-n K )D^\jK)d+ T (j)d VT (j)} 



(51) 



for the terms diagonal with respect to ij. Clearly only the latter ones survive 
the action of the projection operator 
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D.1.5 Kinetic exchange and spin polarization 

In this case of an «inert» /-system the denominator for the resolvent eq. (|T4"1) 
simplifies to: 



where the subscripts omission as compared to the denominator definition by 
eq. (H5|) indicates its special form as diagonal with respect to the /-system 
variables. This quantity is thus the relevant estimate for the effective MMCT 
energy as required in works[39] EH 1121 SSI ED] in the general case. One can 
easily check that the principal source of the renormalization of the MMCT 
energy is the electron-hole interaction apparently omitted there. [39l |4"0| 
H2J H3j [60] Including the dependence of the ionisation potentials and 
electron affinities A(j) of the cZ-shells on the spins of the states produced 
and after assuming the smallness of the intrashell exchange K^d yields the 
final formulae |39[ |4"0"| |4"2"| |4"5] used for semiquantitative analysis. We do not 
expand either LMCT, MLCT, or MMCT energy denominators and use them 
for calculation in their precise form. 

Further derivation evolves as follows: combining the energy denominator 
eq. (|52|) with the projection operators relevant to the considered states in 
the ImH subspace results in the explicit form of the resolvent: 



The hermitean conjugate of eq. (I48|) multiplies the above expression on the 




(52) 




42 



left and after integrating over the /-system variables yields: 



[$ L |vflg)£jw(«)0flgJ?r| $ L > ~EEE 

l^j jJLU pjv' 

E/Wi)/W^ (1+ ^KJ fe^WA-ACf)^ 1 -^^) (i - t»a)J ®ft(&) 

E <v(i)^v(i)OVOid+ (i)^(j) (53) 

tt' 

^M)Mi)D^-\iK) (1 - n„)J ^E^A(0^A(i)^ (1+) 0'A)n A j ©ft (c) 

The terms of the types (a) - (d) represent the contributions[43j used there 
in a simplified form and respectively describe the effects of: (a) transfer 
of holes between the d-shells of the i-th and j-th TMIs through the occu- 
pied /-MOs; (6) and (c) «cyclic» process consisting of transferring an elec- 
tron through the empty /-MOs and of the hole through the occupied /-MOs; 
(d) transfer of electrons through the empty /-MOs. While generalizing the 
configurations[43] these expressions have a sign opposite for the «cyclic» pro- 
cesses (b) and (c) which in their turn took care about the phase relations 
between the MOs in the CN~ anions and more sophisticated bridges. These 
expressions count for the contribution of every pair of the /-MOs into the 
transfer process and also take care about the precised form of the energy 
(denominators) for the MMCT states. [39 j HOI H2J H3] Corresponding config- 
urations are depicted in Figs. [T0 
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D.1.6 Spin-correlated double ionization 

Further contributions of the fourth order in H r come from the ionizing 
components of the effective operator H rr . These are the averages of the 

form: PHrr ± ^QR(E)QHrr T ^ P Corresponding configurations are 

depicted in Figs. [51 [6j The configurations of that type (specifically with 
the two electrons extracted from an /-MO, but not those with two electrons 
added to one of them) have been used,[42j but not systematically. We as 
usually depart from the action of the perturbation operator on the ligand 

(2 + ) 

ground state For Hrr it reads: 

H? r +) \*l) = E E ^PU3)^ 1+ \ 3 K)d^dtA3)a^a KT |$ L ) . 

fiiXcr i/Jkt 

Multiplying this by the resolvent which in the ImQ( 2+ ^ subspace contains as 
multipliers the projection operators to doubly ionized ligand states 
axa'dn'T' a t'T ,a l'a'- One can expect that the singlet and triplet 

ionized states contribute separately. We, however, following|42| concen- 
trate first of all on the states where each /-MO is doubly ionized in either 
sense. That is we take into account only the terms with A' = k'. In this 
case t' = —a' holds and inserting the corresponding energy denominator 
D^ 2+ \\'k' — > ij) = D( 2+ \k'k' — > ij) which is given by 

D (2+) (\k ->■ ij) = I x -A(i)-g x (i) 
+ I K -A{j)-g K {j) 

- 9k (i) - 9x (j) + 9n\ + 9ij (54) 
D {2 -\i 3 ^\k) = I(j)-A x -g x {i) 

+ I{j)-A K -g K {j) 

- 9k 00 - 9x (j) + 9k\ + 9ij 

we arrive to 

E E EE^'^^v(j'O^A«^(j)^ (2+) (Av^ij) 

fi'i'X'a' v'j'k't' /j,iX<T vjkt 

(55) 

which contains the averages of the form (<&l \ a t'T ,a X'a' a ^ a Kr\ $l) over the 
one-determinant /-system ground state which are equal to 
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(3\\'3kk' $aa' $tt' — ^kX'^Xk'^tct'Sctt') n n n x- Inserting this latter value of the av- 
erages and reducing the summation we arrive to the generalization of the 
expressions [42J suitable for programming: 

E E E [/W*)/Wj)f D (2+) (kk ij)DM(jn) [DM(jK) + DW{iK)] (e) 

fii uj K 

(56) 

which contributes to the same exchange matrix elements as the Weihe-Giidel 
cases (i) and (ii) (see Table H]). 

D.2 MMCT contributions to the effective spin operator 

Finally we are in a position to get the estimates of the contributions to the 
effective exchange parameters of the spin Hamiltonian which is in a way 
an approximation to the exact eq. fl53|) . We do this in a line with the 
prescriptions of Ref. (ID] which are based on the fact that if it goes about 
a binuclear i. e. the simplest possible PTMC (a dimer in the terminology of 
Refs. |39 | |4 ^ |4^ | |4^]) each of the terms in eqs. (153|) . (155]) yields a contribution 
to the energy which is proportional to S(S + 1) where S is the total spin 
of the dimer. These contributions when summed up produce the sought 
coefficient at the product SiS-j which is the effective exchange parameter. 
The situation in a PTMC with more than two TMI's is more complicated; 
we consider it elsewhere. For the dimer the derivation evolves as follows: 
each of the terms in eq. (153|) describes a process in which an electron is 
taken from one of the rf-orbitals of one of the TMI's and put to one of the 
(i-orbitals in another TMI. These orbitals are run over by the subscripts fi, v. 
It is only possible in high-symmetry systems when an electron taken from 
one orbital finally lands in another orbital degenerate with the original one. 
However for PTMC (unless it goes about the cuts from highly symmetrical 
solids) one cannot count on the representaion with the dimension higher 
than two. We, however, assume that the symmetry is low enough so that 
no spatially degenerate states can appear. With this precaution we rewrite 
the relevant part of the effective operator eq. fl53l) . The expression eq. ( |53|) 
represents in fact an effective interaction between electrons in the (i-shells 
of the PTMC. The terms off-diagonal with respect to orbital indices fifi' or 
vv' generate correlated transitions between the different crystal field states 
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in the respective (f-shells. However, the magnetic interactions which are 
our main goal and thus the crystal field excitations are not welcome. In 
fact these terms lead to various anisotropic contributions, but we for the 
time being concentrate on the isotropic ones. In order to avoid excitations 
other than magnetic in the ground state manifold one has to set (fi'i) = 
{lii) ; (/i'j) = (vj) ; (v'j) = (vj) ; (v'i) = (/-m) in the above expression and to 
change (fii) -h- (vj) in the terms of the type (b) and (d). 

<$ L |qjflg)m(«)jQflg)*p| $ L > « E 

"£p» K (i)(3Uj)D( 1+ \jK)n K j (^^ x (i)^(j)D^(jX)n^j (a) 

J2KMdurij)OlQid: T (j)d, T (i) 

J2P» K (j)^D^(iK)n)j ^^(i^^D^UX) (1 -n A )J ©^(6) 

(57) 

^<,a)^(i)o;oid+(i)^ r (j) 

- (j2^(mn(i)D^\iK) (1 - n„)j ^/3, A (z)/3, A (j) J D (1+) (jA)n A j ©ft (c) 
J2dt'^Aj)Oi& + dt T (j)d, T (i) 

tt' 

+ (J2M^Mj)D^-\jk) (l-n K )) fefaWhxW^VX) (l~n x )) ©ft(d) 

Y, d Uj)d,T'm\oP_d; T (i)d UT (j) 

tt' 

Further classification is based on the fact |4U [ W2 \ |4"5] that averages of the terms 
entering the above expression depend on the total spin state of PTMC and 
on the way they are composed of the states with definite spins of the involved 
i-th and j-th ^-shells and through this from the occupancy relations between 
the involved orbitals. The required calculations based on the work|76j had 
been performed. |4U| Specifically, each of the terms in eq. (15 7p gives rise to 
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an electron redisribution process (MMCT) starting and ending in a state 
of the total spin S with i-th and j-th e?-shells occurring in the states with 
the spins S(i) and S(j). In the intermediate MMCT state belonging to the 
configuration subspace ImO i f O'L -< Im0 the spins of the involved ci-shells 
can have only one of four combinations of the spin values S(i) ± |; S(j) ± | 
which however may realize or not depending on the availability of electrons 
to be involved in the specific MMCT process or of free space in the d-shells 
at hand. The results are given in Table H] which is derived [4"0] under an 
additional assumption that the total spin of an individual <i-shell always 
goes down whenever a half-filled orbital is involved in the MMCT process at 
this particular shell. 

Otherwise inspecting eq. (j57p shows that for each pair of the d-states 
\xv located respectively in the i-th and j-th TMIs the electron/hole transfer 
through the empty/occupied ligand MOs involves formally the same product 
of the resonance integrals: 

which depending on the type of the process (a) - (d) is complemented by a 
different combination of /-MO occupancy numbers and energy denominators: 

D^\jK)D ( i l J D^\jX)xn K n x (a) 
-D^{j^%D^-\iX)xn K (l -n A ) (b) (58) 
-D^iiK^f^D^UXjx (l-n K )n x (c) 

D^MQ^D^Wx (1 -n x ) (l-n K )(d) 

Either of the MMCT processes (a) - (d) as classified according to the ligand 
MOs involved contributes to the cases (i) - (iv) classified according to the 
type of the MMCT process involved. This allows to establish the correspon- 
dence between the general form of the contributions to the effective exchange 
parameters and their specific forms and cases as given in Table |SJ 

It is remarkable that the overall effect of these contributions cannot be 
represented as one of the effective hopping of electrons between the d-shells. 
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Table 4: Spin dependent factors for the numbers u, L of unpaired electrons in the z-th <i-shell. 



Weihe-Giidel 

cases 


Change of occupancies 

M M 
in the MMCT process 


Unpairity 
factors 
U(ui,Uj) 


Spin dependent 
multipliers 


i 


[1/2] =► [1/2] 


2 


+ 1 5(0 5(0 - 1/2; S(j) -»> S(j) - 1/2 


ii 


[1/2] =► [0] 


2 

Ui( Uj +l) 




[ (+ i) S (i) 5(<) - 1/2; 5(j) " 1/2 
I (_i) 5(0 - 1/2; 5(j) -> S(j) + 1/2 


iii 


[1] =► [1/2] 


2 

(Ui + l)Uj 




[ (+1) 5 (0 ^ 5( - 1/2; S(j) 5(j) - I/ 2 
^ (_1) 5(0 5(0 + 1/2; S(j) 5(j) - 1/2 


iv 


[1] =► [o] 


2 

(Ui + l)(Uj+l) 




[ (+ i) 5(0 _> 5(0 ± 1/2; 5(j) 5(j) ± 1/2 
[ (_1) 5(0 5(0 ± 1/2; 5(j) -> 5(j) =f 1/2 



Table 5: Correspondence between terms [4"2"j |4"5] and eqs. (15 6p . fjoTj) . (158)) ■ 



Factor in Refs. g2l S3] 


Factor in this work 




nxPf.xitj or n K l3^ K {i) 


T/* 


\\ - nx) P,ua( z J or I 1 _ n KJ PvkK\ l ) 


A A 


D^ 1+ \iK)n K 


A* A 


D^(i«)(l-n A ) 


Uab 




v b 


n\/3 u \(j) or n K f3 flK (j) 


v b * 


(1 - 7i A ) /0„aO') or (1 - n„) /VcO') 






A* B 


D^tfA) (l-n A ) 


U BA 









D.3 Correlated double ionization contributions to the 
effective spin operator 

These terms stem from the averages P H^ 2 r +) R^ 2+ \E)Q { - 2 ^ H? r +) P $ i 

and are given by eqs. (155|) . (156)) . The terms within the (2— )-segment are 
processed analogously and result in 
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Figure 1: Sequence of configurations relevant to the matrix elements of type 
(a) as contributing to the effective exchange interaction in the case (i). The 
configuration in the upper left corner belongs to the model subspace ImP^J. 
The configurations in the upper right and lower left corners belong to the 
outer subspace Im<5 < - 1+ - ) - that formed from the model space configurations 
by the LMCT states/processes, that in the lower right corner belongs to the 
outer subspace ImQ - that formed from the model space configurations by the 
MMCT states/processes. Transitions between the subspaces are indicated by 
the arrows, whereas the specific electron transfers by the arc arrows. The 
contributions to the exchange parameter are obtained as products of the 
quantities indicated in the Figure. The case (i) is characterized by the choice 
of the values of the energy denominators: namely those with the lowering 
spin in the both z'-th and j-th <i-shells. 
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Figure 2: Sequence of configurations relevant to the matrix elements of type 
(b) as contributing to the effective exchange interaction in the case (i). The 
configurations in the upper right and lower left corners belong to the outer 
subspaces hxxQ^' - that formed from the model space configurations by the 
LMCT states/processes. Other notation is as in Fig. [TJ 



52 



W K -t-^^O) 
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II 



7C /C 



— (- »/> at) (M) -41- P <= — -H-> 



Figure 3: Sequence of configurations relevant to the matrix elements of type 
(c) as contributing to the effective exchange interaction in the case (i). The 
configurations in the upper right and lower left corners belong to the outer 
subspaces ImQ^ lT ^ - that formed from the model space configurations by the 
LMCT states/processes. Other notation is as in Fig. [TJ 
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Figure 4: Sequence of configurations relevant to the matrix elements of type 
(d) as contributing to the effective exchange interaction in the case (i). The 
configurations in the upper right and lower left corners belong to the outer 
subspaces ImQ^ 1- - 1 - that formed from the model space configurations by the 
LMCT states/processes. Other notation is as in Fig. [TJ 
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Figure 5: Sequence of configurations relevant to the matrix elements of type 
(e) as contributing to the effective exchange interaction in the case (i). The 
configurations in the upper right and lower left corners belong to the outer 
subspace lmQ^ 1+ ^ - that formed from the model space configurations by the 
LMCT states/processes, that in the lower right corner belongs to the outer 
subspace lmQ^ 2+ ^ - that formed from the model space configurations by the 
(LM) 2 CT states/processes. Other notation is as in Fig. [TJ 
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Figure 6: Sequence of configurations relevant to the matrix elements of type 
(f) as contributing to the effective exchange interaction in the case (i). The 
configurations in the upper right and lower left corners belong to the outer 
subspace ImQ^ 1 ^ - that formed from the model space configurations by the 
MLCT states/processes, that in the lower right corner belongs to the outer 
subspace ImQ^ 2 ") - that formed from the model space configurations by the 
(ML) 2 CT states/processes. Other notation is as in Fig. [TJ 
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